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Abstract 



It is shown some exact solutions in the Brans-Dickc (BD) theory for a 
Bianchi V metric having the property of inflationary expansion, graceful exit, 
and asymptotic evolution to a Friedmann-Robertson- Walker (FRW) open 
model. It is remarkable that an inflationary behaviour can occur, even without 
a cosmological potential or constant. However, the horizon and flatness prob- 
lems cannot be solve within the standard BD theory because the inflationary 
period is severely restricted by the value of the BD parameter u. 
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1 INTRODUCTION 

It is well known that the Universe is homogeneous and isotropic on very large scales. 
Supporting this assertion are different measurements as the isotropy measured in 
the Cosmic Microwave Background Radiation (CMBR) by the Cosmic Background 
Explorer (COBE) satellite JTJ; isotropy in x-ray backgrounds (e.g., quasars at high 
red shift); and in number counts in faint radio sources, cf. Ref. ||. All these 
measurements give evidence for the cosmological principle, which consistently should 
be valid on very large scales of our Universe. However, on smaller scales there are 
a variety of inhomogeneities and anisotropies that have been produced by grow of 
matter perturbations, at least since the time of last scattering until now. Accor- 
dingly, many authors have tried to model the early Universe with non-FRW metrics 
allowing general initial conditions to explain, after its evolution, small and large 
scale properties of our Universe. Thus, the chaotic cosmology programme || tried, 
unsuccessfully, to isotropize anisotropic cosmologies in order to understand, within 
general footings, the present large scale isotropy. Collins and Hawking || proved, 
within General Relativity (GR), that if the dominant energy condition and positive 
presssure criterion are satisfied, the universe can approach isotropy only if it is one of 
the types I, V, VIIo, and Vllh- However, at finite times other anisotropic models can 
be bounded from above || and even they are consistent with COBE measurements 
H, that is, anisotropic cosmologies do not need inflationary scenarios to explain the 
observed ~ 10~ 5 in the CMBR. However, a causal origin of the perturbations 
is not provided by anisotropic cosmologies alone, and inflationary cosmologies does 
provide such an explanation, up to some fine tuning. This is one reason why inflation 
is still the most interesting proposal to the standard model of cosmology. 

The question of the entry to an inflationary phase, now known as the cosmic 
no hair conjecture 0, as well as the question of its exit (graceful exit problem) are 
two key aspects of inflation that have been intensively discussed. No hair theorems 
have been proved in the context of GR within anisotropic and some inhomogeneous 
cosmologies ||. The entry to inflation has been also discussed in some other con- 
texts, including tilt matter models ||, and some anisotropic models in scalar tensor 
theories [|Hj] . In general, one has realized that for some set of initial conditions there 
exists an inflationary attractor for the evolution. However, the exit of inflation is 
not necessarily guaranteed when a cosmological constant (A) is present. Instead 
of A one has been considering a cosmological function V(<f)) that plays the role of 
a constant during some time r = NH~ l , where N is the number of e-foldings of 
exponential expansion and H ^const. the Hubble parameter. Adding a properly 
chosen constantQto V(<p) one gets after inflation that V(<ftf) ~ 0, then the Universe 
1 So fine tuned to avoid having a cosmological constant nowadays (n) bigger than 8irGp n = 
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preheats [|TT], reheats, and experiences a metamorphosis to a FRW dynamics filled 
with a particle content characterized by the decay of the 0— particles. The intro- 
duction of a constant in V(4>), to allow inflation to occur and to graceful exit, is 
well motived but it is put by hand. In this sense, it would be desirable to have a 
model in which both features, the entry and exit of inflation, arise naturally, as a 
consequence of the exact, integrating dynamics. In GR, however, one needs to in- 
troduce a cosmological constant (or function) to generate inflation. In scalar-tensor 
theories this term appears naturally, but to obtain successful inflation it implies the 
fine tuning of parameters that not always are in accord with the particle physics 
models behind. 

Recently, effort has been put on the analysis of FRW models in the Brans-Dicke 



(BD) theory [0, which is simplest gravity theory beyond GR. In Refs. |T3|, [14| a 
qualitative analysis of the general evolution and asymptotic behaviour is studied, as 
well as the conditions for inflation even without a cosmological function or constant, 



see also Ref . |15j . Motived by these analyses we study in the present work Bianchi V 
anisotropic solutions in the BD theory, because it is one of the simplest step further 
in complexity and permits the analysis of the anisotropic properties of models; in 
particular the isotropization of some solutions has been studied in Refs. [[LJ], [T7|j . 
We found that Bianchi V models can have the properties of inflation (without a 
potential or a cosmological constant), graceful exit of the inflationary stage, and 
evolution to an effective FRW open Universe, like the one we observe nowadays. 
The range of values for the coupling constant uj is, however, very different from 
what one would desire, i.e. to > 500 



In the next section we show the analytic solutions and, thereafter, in section [| 
we analyse their behaviour. In section |4] some comments are made on fluctuations 
of the BD field in these models. Finally, section |5] is devoted to conclusions. 

2 BRANS-DICKE BIANCHI V SOLUTIONS 



The BD theory JL2] has the following Einstein-Hilbert action, with signature (+,-,- 
,-): 



(1) 



c 



>R+ + 16?rL 



M 



16tt 

where R is the Ricci scalar, the scalar field, the symbol | partial derivative, u the 
coupling constant of the theory, and g the determinant of the metric tensor. After 



ml ~ 10" 83 GeV 2 . 
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varying this equation one derives the BD field equations 



R 



Rg 



-T, 



uJ 



Wv - IIA 9iw 



(2) 



and 



A 



8tt 



where the symbol 1 1 stands for the covariant derivative and T is the trace of the 
energy- momentum tensor, T^. The continuity equation (energy- momentum con- 
servation law) reads 

V 11* = • w 

We consider homogeneous, anisotropic Bianchi-type models to study the dy- 
namics, and particularly a Bianchi type V spacetime symmetry in a synchronous 
coordinate frame given by the line element [19]: 



ds = dt — a l dx — a 2 e x dy 



a 3 e 



-2x 



dz 2 



(5) 



Following, we shall use scaled variables to find and analyse the solutions in a simple 
way: the scaled field ip = <pa 3<yl ~ v \ a new cosmic time parameter dr\ = a~ 3u dt, ()' = 
J^, the 'volume' a 3 = aia 2 a 3 , and the Hubble parameters Hi = a//aj corresponding 
to the scale factors a* = aiirf) for i = 1,2,3. We assume comoving coordinates 
and a barotropic equation of state for the perfect fluid present (T^), p = up, with 
v a constant. Using these definitions and the above given metric, one obtains the 
cosmological equations from Eqs. (0) and (0): 



8tt 
3 + 2a; 



[1 + (1 - v)u\ pa 3(1+l/) for i = l,2,3. (6) 



if;' 



H X H 2 + H X H 3 + H 2 H 3 + [1 + (1 - u)u] (#1 + H 2 + H 3 ) 



1 -v)[l + u{l - v)/2](Ht + H 2 + H 3 

,3(1+*) 



8^ 



pa" 



(7) 



i[>" + (v-l)a 6v C v il) 



57T 



3 + 2cj 



[2(2 - 3i/) + 3(1 - z/) 2 w] pa 



3(1 + !/) 



(8) 
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and 

H 2 + H 3 = 2H ± , (9) 

where Cy = T.iCiv with the partial the curvature terms equal to each other, i.e., 
Civ = 2/af for i = 1,2,3. Eq. (pj) implies that a>2 and a 3 are inverse proportional 
functions, 0,2(13 = af; note that 3H = Hi + H2 + H3 — 3Hi, H being the mean 
Hubble parameter. 

Additionally, the continuity equation yields: 

pa 3(1+l/) = const. = M v (10) 

My being a dimensional constant depending on the fluid present. The vacuum case 
is attained when M v = 0. 

The system of ordinary differential equations, Eqs. (§-|D, can be decoupled to 
obtain an equation for ip alone, by means of the relation a,\ = a, which can be once 
integrated to get: 

^ v - 37rby^ 2 _ 2 3(i~-u) K(1 _ 3u)v + Vo] ^ + p + C 1 - + 3 "VH 

Of 1 _ liA m 

3 (1 _ ^ [2(2 - 31/) + 3(1 - v)M [^(1 - 31/) 77 2 + % 17 + ^]m„ = , 

where r) , 771, and -K" are integration constants; m„ = fr^j. Solutions to Eq. (]TTf) 
are to be used to get the mean Hubble parameter through the following equation 



Hi + H 2 + 



1 - z/ 



V>' (1 - 3v)m v r} + r? 



(12) 



obtained by using Eqs. (0) and (H). Additionally, each Hubble rate is written as 
follows (similar to the Bianchi type I model Ref . fl2(J ) : 

^ = 3(^1 + ^ + ^) + ^ = 3^3^ , (13) 

where the h^s are constants that determine the anisotropic character of the solutions. 
If hi = for % = 1, 2, 3 simultaneously, no anisotropy is present. Furthermore, the 
Bianchi V model obeys the condition 

hi + h 2 + h 3 = 0, with h t = 0, (14) 

the last relation stems from the fact that Hi = H, see Eq. (||). 

In order to analyse the anisotropic character of the solutions, we have constructed 
the following 'constraint' equation for the anisotropic shear, a: 

<t(ti)= -(Hi-H2) 2 -(H2-H 3 ) 2 -(H 3 -Hi) 2 = 
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>"\ 1 fip'V (l-3j/)/(l-3i/)m„i)|)) \ /f 



[2 - 3v + |cu(l - z/) 2 ] / (l -3^)m^ + ?7 \ 2 3[2 + iu(l - z/)(l + 3z/)]m 
+ (l-^) 2 V i> ) + 2(1 - v)i) ^ } 

a = is a necessary condition to obtain a FRW cosmology since it implies Hi = 
H2 = H3, cf. Ref. [ 2"I| , |TE| |. If the sum of the squared differences of the Hubble 



expansion rates tends to zero, it would mean that the anisotropic scale factors tend 
to a single function of time which is, certainly, the scale factor of the open Friedmann 
model. 

The anisotropic shear becomes, using Eq. fll4]) , 

= -^r . (16) 



or the dimensionless shear parameter | 



a 54(1 - vfh 
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H 2 [ip' — (1 — ?>u)m u rj — r] ] 2 ^ 

These equations admit solutions such that a — > 0, a/H 2 — > as r\ — > 00 (or 
t — > 00), that is, one has time asymptotic isotropization solutions, similar to the 
solutions found for the Bianchi V model in GR, see Ref. ||. In fact, one does not 
need to impose an asymptotic, infinity condition, but just that r] ^ 77*, where 77* is 
yet some arbitrary value to warrant that a (a/H 2 ) can be bounded from above. 

The above equations are valid for any value u of the equation of state, but the 
solutions are different for v ^ 1/3 and v = 1/3. Therefore, we treat both cases 
separately. 

2.1 Solutions with v 7^ | 

The explicit solutions are following given. A particular solution of Eq. ( |TTD is |16|| : 



if> = Avtf + Byri + Cv, (18) 



with the constants 



3u 



|2 



A " ~ (1 + 3 „) m - ' 

^ = -Klrf)"- (19) 

being the same as for the isotropic, k = ±1 cases|16[], but now 

m u [l + 3v)C v - 18u + u r 1 + 3u) 2 ■ (20) 
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The constants B v and C v , being proportional to the h^s, encode information on the 
anisotropic shear, that is, on the nature of the anisotropic character of this Bianchi 
type model. The constants rji and K of Eq. flTT|) are determined through Eq. (|2"0| ) 
but they have no further significance. 

The Hubble expansion rates, given through Eqs. ([13|, [Tj|, [19]), are 

1 (1 - 3u)m u r] + rjo — (1 + Su)hi 



Hliv) (1 + 3//) {A v rf + B v r) + C v ) 
which can be directly integrated to get the scale factors: 

-2(1 + 3u) W^l r 



for i = 1,2,3 (21) 



a i 



[2 + (l-i/)(l + 3i/)w]m„ 



A v tf + B vV + C v 



2(l-3i/) 



or 



a 2 = ai exp 



a 2 = a\ exp 



-2ho 



A 

2h 2 



■arctanh 



:arctan 



-2(1 — 3^)[(1 - 3v)m u ri + V o}' 
-2(1 - 3i/)[(l -3z/)m^ + r7o] N 



and according to Eq. @ one gets 



;i + 3z/)v^A 



A > 



A < 



(22) 
(23) 

(24) 



a 3 = — 
a 2 



where the discriminant is 

A = B 2 , - AAC 



-8(1 -3v) 



18u + 1 + 3u ) 2 uj 



(25) 



(26) 



being proportional to the shear, see Eq. (p!q). Independent of the value A might 
have, hi + h 2 + h s = is always true. The type V model with A ^ fulfills that 
h 2 = —h 3 with hi = to have truly anisotropic solutions. But for A = 0, 



O3 — Oj 2 — &i 



(27) 



since one has that hi — h 2 — /13 = 0. In this case, the scale factors are equal to each 
other, given by Eq. (p2|), up to a constant that can be scaled away. 
The BD field is 



a 



3(l-«0 



[2 + (l-i/)(l + 3i/V]m„ 
-2(1 + 3z/) 



3(1-./) 

2(l-3i/) . l+3i/ 



(28) 
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The above-presented solutions reduce to the previously known for the isotropic, 



FRW model with an open space (k = — 1), cf. Refs. |f23|, ^4], when A = 0, 
implying that hi — hi — hs — 0, to have no shear, a = 0. The above solutions have 
no analogues in GR since the scalar field <fi cannot be set to a constant. 

It has been conjectured, for the Bianchi V dust [y = 0) model, that it is almost 
impossible to find solutions because of the complexity of the equation for ip, which 
is presented in Ref. |26| as a fourth-order differential equation. Therefore, only 
particular solutions are expected to be found, if any. In the present work it is 
shown an integrated equation for ip which is now of second-order, see Eq. (|TTD , 
and possesses at least one solution for v ^ 1/3 given by Eq. (|i~8D ; other solutions 



are hoped to be found because of the relative simplicity of Eq. (O). The dust 



solution in Ref. p6| is written in other variables, but it is a special case of our above 
solution when one takes B — v — 0. Referring to this, one of the main problems 
of the physical interpretation in the field of exact solutions in cosmology has been 
that results are written in scaled variables and the physics behind them is partially 
hidden. In our case, for the dust model [y = 0) the scaled time is the same as 
the original cosmic time, written in Eq. (|5]). For v ^ effort should be made to 
interpret physical situations. 

Following we present the solutions for the radiation case. 

2.2 Solutions with zy = 1/3 

In the radiation case the curvature terms in Eqs. (|| [7], ||) turn out to be especially 
simple (a 6u Civ — 2, a 6u Cy = 6). Therefore, Eq. (§) reduces to 

167rilf 

V"-4V> 3"^ = °' ( 29 ) 



2 



that can be directly integrated to get 

47rM 

ip = C\t v + c 2 e v — , ci, C2 are integration constants. (30) 

3 



The general solution for the Hubble expansion rates, using Eqs. (|i~2j [13], |3(|), is 
given by 

Hi{n) = ^ for 2 = 1,2,3 (31) 

c x e~ 2r > + c 2 e 2r i 3^ 

where C3 (=. rj /2) is an arbitrary integration constant and the constants hi, account- 
ing for the anisotropy of the models, are determined by: 

h = 0, 



This solution is also valid for Eq. (|ll| ) with v = 1/3. 



9 



h 2 



=Fa/12A - (3 + 2u) c| 



(32) 



where A = (27rM 1/3 /3) 2 — c\c 2 . The constants given in Eq. (|32|) are consistent with 
Eqs. (|1|, ID). For if h x = h 2 = h 3 = 0, it implies that cx = 0. 
From Eq. ([31]) one obtains the scale factors 



a io ip 1/2 exp 



c 3 - K 



cii = a io ^ 1/2 



or 



exp 



CI; 



. 2VA 

fri - c 3 



arctanh 



c 2 e 2 " - 27TM/3) 



2 V /Z A 

= a io ^ 1/2 



arctan 



c 2 e^-2nM 1/3 /3) 



exp 



c 3 - h 



, A = 



A>0, (33) 
A < , (34) 
(35) 



_2c 2 e 2 '? -4vrM 1/3 /3_ 

where «j are integration constants proportional to some initial size of the Universe. 
Note that Eq. fl2~5|) is valid again. 



The BD field (0 = ^/a 2 ) is 



{a lo a 2o a 3o ) 2/3 exp 



(ai a 2o a 3o ) 2/3 exp 



^3 



A 

c 3 



arctanh 



c 2 e 2 ^-27rM 1/3 /3) 



x^A 



arctan 



:c 2 e 2 ^-27rM 1/3 /3) 



A > , (36) 

, A<0, 

(37) 



or 



(ai a 2o a 3o ) 2/3 exp 



-c 3 



A = 



[oe*! _ 27rM 1/3 /3j ' " (38) 

From Eq. (|32]) one observes that if A < 0, then 3 + 2u < in order to have real 
solutions (h 2 , h 3 being real numbers) for the scale factors Eq. (|34|) or Eq. (pop. 

Eqs. (|33| - p8D , together with Eq. (0), show the most general solution for the 
Bianchi model V containing a radiation field or a fluid of relativistic particles. For 
some particular values of our constants ci, c 2 , c 3 , and hi, they were firstly presented 
in Ref. PE| . The solution for « = 1 is the most general solution for the FRW open 
Universe, cf. [27]. In the particular case that c\ = —c 2 and hi = one gets the open 
FRW solution given in Ref. p8| , where a misprint is found: (3 — *> —(3 in Eq. (76) of 
that paper. 

In the particular case that C\ = c 2 = one obtains the isotropic solution with <fi < 
0. In this case, one has that H(rj) =const., ip(t]) =const. and since a(r])dri = dt, one 
has that Hi(t) = H 2 {t) = H 3 (t) = 1/t with a\ = ai D t, a 2 = a 2o t, a 3 = a 3o t, getting 
that (t> = — „ji L , a sign of antigravity?, a particular solution found elsewhere |29 
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For if C3 = 0, the solution becomes the same as in the GR theory with =const= 
Mp h first obtained in Ref. f50fl . Further, if hi = C3 = 0, one obtains the k — — 1 



isotropic solutions of GR discussed in the late 60's 



3 SOLUTIONS ANALYSIS 

The solutions presented in the preceding section show various behaviours depending 
upon the value of the BD parameter (u), the type of fluid present (u), the amount 
of matter (m u ), and the anisotropic degree (hi) of the models; let us consider each of 



them. In the BD theory measurements imply that uo > 500, cf. Ref. [fL8|| . However 



in general scalar tensor theories the coupling constant can be much less than one 



3^, |33|. If such a theory is to be valid during some epoch in the early Universe, 
then it can be mimicked by an effective BD theory, but with its proper value of uo. 
Therefore, it is interesting to consider both great and small cu-values. 

In the preceding section we have presented two types of solutions for models 
with v ^ 1/3 and v = 1/3. Following, we present the analysis of these solutions 
separately. 

3.1 The case v = 

Within the most interesting barotropic equations of state (p = up) with v ^ 1/3 are 
the cases of dust (v = 0), stiff fluid or supernuclear density (v = 1) |34|], particle 



creation (v = 2/3) [|35|j , and vacuum energy (v = —1). For concreteness we shall 
consider in this subsection a gas of nonrelativistic massive bosons that could have 
played a dominant role (in the stress energy tensor) for the dynamics of the very 
early universe. This gas behaves as dust with an equation of state determined 
by v = |36|]. Accordingly, the total mass of the gas is given by the quantity 
M = ^m = p al 

The values of hi are related to the degree of anisotropy present in the model. 
As pointed out already, if hi = then A = implying that (X\ — Qj2 — O3, to 
have the open FRW model. Anisotropic solutions are found for the case that the 
above constants fulfill that A ^ 0. For if A > with v = 0, then u < 0, cf. Eq. 



6f). In this case the solution is given by Eqs. (|22|, 23, pa) implying that it is valid 



only in the intervalPl V 2 I UJ IM vo y/ 2/w \h 2 \ no ^ <jY ien f or fais solution to 

J T-l m 1 m 

evolve long times one should demand that 3> that is, this would be 

favored by large anisotropics (h%), small masses, or small negative a;- values. This 
solution is for hi > deflationary in the scale factor a>2 and inflationary in a^, for 

3 Note that for v = one has that r\ = const, t. 
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h,2 < the roles of a 2 and 03 are inverted. However, there is no exit of inflation 
because of the continuously increasing arctanh function in Eq. (|23|). Therefore, this 
behaviour can be only useful if after some inflation time some other physical source 
(e.g. another scalar field) begins to dominate the dynamics to stop the otherwise 
eternal inflationary behaviour. 

The most interesting case seems to be for A < 0, which is actually the case for 
u > 0. The solution is given by Eqs. (|22| , p4| , [2~5| ) and is valid in the whole interval! 
— 00 < 77 < 00 without restrictions. For h 2 > this solution is initially inflationary 
in ci2 and deflationary in a^; again their roles can be inverted by choosing h 2 < 0. 
In contrast with the A > solution, in the present case the solution has exit, since 
the arctan function in Eq. (|24]) tends asymptotically to ~. Thus, during some time 
interval this function behaves inflationary, after which a FRW behaviour follows. 

The strong exponential expansion is possible, though no cosmological constant 
or function is present; this is because the non-minimal coupling in Eq. ([!]) implies 
new kinetic terms to the dynamical equations in comparison with GR ||15|| . In 
order to achieve a successful inflationary scenario one should demand that the three 
scale factors, or its volume, inflate simultaneously. This is apparently not possible 
from the form of the solution for ai, Eq. (|22|), which is a power-law standard 
expansion type. An inflationary stage may occur, however, when the denominator 
of Eq. (^) is effectively a constant (ip ~ const.), while the numerator is a linear 
function, so Hi ~ 77. This behaviour can be attained only during some time interval, 
since ip, given by Eq. flTBp, is a quadratic function and is eventually numerically 
larger than the numerator to give rise to Hi ~ 1/rj. To get enough e-foldings 
(^- = e N ) of inflation to solve the horizon and flatness problems, one can try to 
adjust the model parameters to, hi, and m . Let us impose the necessary conditions 
for it. On the one hand, ip ~ const, from rj = until the time 77*, with 77* = 

^— (vn + — 1 — — • On the other hand, to last sufficient amount of e-foldings of 

mo V ^ J m 

inflation^ (N ~ 68) the accelerated expansion must run at least until the time 77/, 
with 77/ = — [(?7o - hi) 2 + 2N(r]% + ^)j 1/2 - Thus, the inflationary epoch is 



valid in the interval < 77 < 77/, provided that 77/ < 77*. Unfortunately, the latter 

2/j2 

inequality is valid for very restricted values of a;, in fact for uj < Otherwise, 

for arbitrary uj the number of e-foldings is at most N — 1/2 only. For instance, 
one may try to fit the parameter m Q (if m is augmented, N grows at a fix time) 
to get the desired number of e-foldings, but then the above inequality is no more 
valid. Analogously, whereas arbitrary large values of — can make arbitrary large 



the inflation time ||37|| , the time 77/ grows in the same proportion, avoiding to have 

4 For initial conditions with lower energy scales than Planckian, the number of e-folds can be 
smaller than 68 0|. In our models we pursue to identify the asymptotic value of the scalar field 
with the Newtonian constant G. Therefore, the energy scale G 1 / 2 appears. 
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T]f < 77*. In BD FRW cosmologies a similar result is found, where only negative 
values for to can accomplish an accelerated, successful expansion |I5fl: A solution to 



the flatness problem is not possible because the current value to = 500 [H| resctricts 
the solution types, not allowing an accelerated expansion, and not allowing to have 
nowadays (n) a value of 0" 1 = IQttG, G being the Newtonian constant []38f| . 

In figures [l|-|3] analytic solutions are plotted for some cosmological parameters. 
Figure [I] shows the Hubble parameters and the BD field for v = 0, uj = 500, and 
h 2 = Vo — m o — 1- Initially the Hubble parameter H 2 grows almost linearly 
until 7] = 0.42 and, therefore, the solution is inflationary in this direction with 
a 2 ~ e Kir?+K2r?2 , where k±, k 2 are constants. Afterwards H 2 continues growing at 
a lower rate until it stops at rj = 1.0 . After this time this Hubble parameter is 



dominated by its denominator, cf. Eq. ([13]), and the solution becomes H 2 ~ -. 



Later on, after rj = 1.6 all three Hi are dominated by its denominator and the 
solution is of FRW type. 

Figure ^| shows the solution for uj = 500, v = 0, r] = m = 1, and h 2 = 10. 
In this case, the initial anisotropy parameter was augmented, provoking an initial 
contraction of a 2 (H 2o < 0), after which an expansion follows. If one leaves above 
parameters fixed, but one augments the value of m , a 2 goes from a contraction to 
an inflation period. 

Finally, figure |3| is similar to figure [TJ, but in figure [3] a small u was chosen to 
mimic a model of induced gravity valid when its potential plays no role, see Ref. 
[ P7| . In figure |3] we observe an initially inflationary behaviour for the three scale 
factors, its exit, and its evolution to the open FRW. Again the same effect is attained 
by augmenting m , the three scale factors inflate but with a restricted number of N 
e-foldings, as mentioned above. 

Note that during the inflationary era the solutions are still anisotropic; they tend 
to isotropize once Hi ~ I/77, see figures [l], H, anc ^ B 

For small cu's, like in some induced gravity modelsf] [^|, the three Hubble param- 
eters grow, as can be observed from figure |3|. For uj 1, cf. Ref. |33|, the models 



isotropize during the inflationary era |37|]; later on, they evolve asymptotically to 
the open FRW model. 



(39) 



The dimensionless shear parameter, Eq. (|I7D, turns out to be 

a 6ho 



H 2 [m v rj -rio] 2 

which for the A < solution, one observes an asymptotic, isotropic behaviour 



5 Our solutions are applicable in the induced gravity theory during a time interval, when its 
potential is not significant for the dynamics p7[. 
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Figure 1: 
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Figure 2: 
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Figure 3: 
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a/H 2 — > 0, when 77 — > 00, to yield the following solution: 




(40) 

the bigger u> is, the larger (smaller) is the ratio between a\ and a 2 (a 3 ). In our 
graphics we observe only small differences in the Hi (i = 1, 2, 3) for 77 ~ 20 to get, 
in the asymptotic limit, the solutions found in Refs. |23], |24], p5[j . This asymptotic 
solution is different from the open FRW solution in GR, and is described by a fixed 
point in the autonomous phase plane (a/a,<p/4>) analysis for FRW cosmologies in 
the BD theory [0 . 



3.2 The case v = 1/3 

The general solution found for v = 1/3 is valid for a fluid of ultra-relativistic particles 
and/or radiation. The isotropic, open solutions are obtained when hi = 0, and for 
some particular values of our constants they are reported in Refs. |27|, |28j . Solutions 
given by Eqs. (p3|)-(p5|) have a similar structure: they are multiplied by the power- 
law term "ip 1 ^ 2 , as well as by the exponential term. Initially the exponential term 
maybe more significant, and from it one could possibly obtain an inflationary stage 
in the evolution. The time asymptotic behaviour is dominated by the power-law 
term. The details depend on the value of the discriminant (A). For the A = case 
the number of e-foldings of inflation is given by 



2c 2 e 2r >f ^ 2c 2 ^ 2 y° 



where r/j is the time when the power-law solution begins to dominate, i.e., when 
the exponential term tends to one. The number of e-foldings depends finely on the 
value of the constants c 2 , c 3 , and M 1/3 , subject to the following inequalities: (i) c 2 > 
27rM 1/3 /3 > 0, the first inequality is to avoid the point where becomes infinity, and 

the second is to have a«, p > 0. Note that c\ > since C1C2 = (2irM 1/3 /3^ > 0; (ii) 
Ci + c 2 > 47rM 1/3 /3 to have if) , a io > 0; and (Hi) c 3 < to have N > 0. The right 
choice of these parameters gives the desired number of e-folding of inflation (N). 
Unfortunately, the present case (A = 0) is valid for uj < —3/2, otherwise a 2 and a 3 
become complex functions, cf. Eqs. (|32| ) and (|3^). Recall that values for u < —3/2 
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correspond to a negative scalar field density in the conformally rescaled Einstein 
frame (g^ = Gcfig^), see for instance [14]]. However, if one believes that the correct 
physics is in the Jordan frame, negative a/s maybe still interesting. Recall that 
changing frames via conformal transformations imply physically different gravity 
theories since the matter couplings are different, as well as geodesic equations, see 



for instance 39 



Solutions with A > have the same problem as in the v ^ | case: they inflate 
anisotropically during the whole evolution given by the time interval 

'2-kM 



In 



1/2 



(2vrM 1/3 /3 - l) -f\' 2 < V < In (2ttM 1/3 /3 + 1 



1/3 _ Ci 



1/2] V2 



and the graceful exit problem remains because inflation never ends. 

Solutions with A < imply that u < —3/2 in order to have real (instead of com- 
plex) solutions for the scale factors in accordance with Eqs. (32) and (134]) . Figures 



|4| and | show solutions for the different values of the physical parameters. One notes 
that terms containing C2e 2r] dominate the dynamics very rapid, and therefore shall 
account for the major contribution to the dynamics in asymptotic times. Terms 
with cie~ 2v could be important for the very beginning. In figure § we have chosen 
the parameters to be u = —2, h 2 = c\ = c 2 = 1, and m 1/3 = —1, which corresponds 
to M 1/3 = +1/(8-71") to have a positive density. In figure [| we have augmented the 
value of the anisotropic parameter h 2 (=100), and the values of c\ and c 2 (ci > c 2 ) as 
well, to observe the influence of the Cie~ 2r? term. Otherwise, if c 2 > c\ the solution 
becomes the isotropic solution yet from the very beginning. 

The solutions with A < have the chance to achieve a number of e-foldings of 
inflation (^ = e^) given by 



N 



-c 3 



2 V /Z A 



arctan^ _ _ arctan * L _ ^ 



1, ip f 
+-ln-^- 
2 

(42) 



One can choose the integration constants to have the desired number of e-foldings, 
subject to following constrains (the arguments are the same as in the A = case): 
(i) c 2 > 0, M 1/3 > 0. Note that A < implies that cic 2 > (2vrM 1/3 /3) 2 > 0, hence, 
ci > 0; (ii) C\ + C2 > 47rM 1/3 /3; and (in) C3 < 0. Accordingly, figure ^| shows the 
parameter dependence of the number of e-foldings of inflation (N) . Small values of 
C2 together with big hi values are requested to achieve a sufficient number for N. 
The dimensionless shear parameter, Eq. (|17D , turns out to be 

a 6hl 

IP = ~ [- Cl e-^ + c 2 e 2 n - c 3 ] 2 ' (43) 
which for the A < solutions, their asymptotic behaviour, when rj — * 00, implies 
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Figure 5: 
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C2 

8 IB 24 32 4Q ijg 




a I H 2 — > 0, yielding the following isotropic solution: 



Hi(r}) -> 1, ^(77) -> a. 



c 2 e" 



0(77) 



(aio a 2o a 3oJ 



-2/3 



-2/3 



(44) 



Independent of the values all parameters may have, the solution evolves to its 
isotropic behaviour ||29|| , as can be seen in figures |] and |^. The above limit so- 
lution, Eq. (H), tends to the radiation solution in GR for asymptotic times. 



4 FLUCTUATIONS OF THE BD FIELD 

Although we are working in a framewok of a classical theory, one may also consider 
quantum fluctuations of the BD field during inflation. It is well known that isocur- 
vature and curvature fluctuations are generated, and they produce perturbations in 
the density and in the CMBR. Fluctuations of the BD field has been analised 0, fyj] 
mostly in the Einstein framed], since there the new scalar field ((f)) is minimally cou- 
pled and its the kinetic term is canonical, therefore, the role of fluctuations is easy 
to interprete physically, and to import results from what is known in the Einstein 



theory j|2], [40| . For example, if it were the case that we consider an additional effec- 
tive constant potential (false vacuum energy) in the Jordan frame, this would turn 
out to behave like an exponential potential in the Einstein frame, i.e., a la extended 



inflation |HJ. In this case, one gets that Sp/p depends on A 4 ^ 2w ^ due to power 



law inflation, where A is some perturbation scale [^3], |4(| [0]. Then, all machinery 
of reconstructing the inflation potential could be applied J|4]. In our case, one does 



not have a potential that forces H = const., but in our models one has that H ~ t 
during some initial stage. Because of this fact, one expects that the spectrum is 
non-scale invariant too. In our inflationary models, however, if one considers the 
Einstein frame to analise the perturbations, one finds that the transformation fac- 
tors involve quantities with ^3 + 2u (see footnote ||), which in our solutions with 
inflation, it is a complex number. Then, one cannot correctly interprete our infla- 
tionary solutions in the Einstein frame. Furthermore, values of ui < —3/2 generate 
negative scalar-field densities in the conformally rescaled Einstein frame [14| . 



For us > —3/2 we do not have an inflationary behaviour, therefore, the production 
of quantum fluctuactions is out of interest because the horizon problem remains. 
However, if one starts the anisotropic models with Planckian initial conditions one 
can explain || the observed anisotropy of the CMBR measured by COBE JTJj. 



3 The transformation of the metric is — Gcftg^ and the new field 4> — \ itwG l n (G</>) 
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5 CONCLUSIONS 



Our analysis is based on the exact solutions reported in this work and in Ref. fl6 
For a fluid with v ^ 1/3 the particular solution Eq. (|l^) was used, but other 
solutions may be found from Eq. (|TT1), which is considerable simpler than other 
equations found in the past [ffijfl . For instance, our time asymptotic solution does 
not tend to the GR one, therefore, the latter must correspond to another solution of 
Eq. For the radiation case (u — 1/3) the most general solution is given. The 

FRW solution is same as the given by our i = 1 case (recall that h\ = 0). The GR 
Bianchi V case is achieved when C3 = 0. 

The above-presented exact solutions have a variety of possible behaviours de- 
pending on the values of the physical parameters u>, z/, m u , and hi- Within some 
parameter range it is possible to encompass a cosmological model that incorporates 
the issues mentioned in the introduction: to have anisotropic initial conditions that 
leads to an inflationary stage, and, as time goes on, to isotropize towards an open 
FRW model. In this way, inflation represents a transient attractor, and the FRW 
behaviour is an asymptotic attractor. Further, these models can explain the ob- 
served anisotropy degree of the CMBR measured by COBE |]J, if the model is to 
be started with Planckian initial conditions ||. 

It is peculiar that these solutions show an inflationary behaviour, even without a 
cosmological constant or function. This is because the non-minimal coupling in Eq. 
([]]) implies new kinetic terms to the dynamical equations in comparison with GR fl5| . 



For the non-radiating [y 7^ 1/3), anisotropic solutions with A ^ the inflationary 
stage takes place when the denominator of Eq. (|T3|) is almost a constant, whereas 
the numerator is linear. This is achieved during some time interval depending on 
the physical parameters u, u, hi, and m v . Especially, large values of m v or — favor 
strong exponential expanding solutions of the type ~ e^i^H-^ 2 ^ ^ut the number 
of e-foldings of inflation cannot be greater than 1/2 for uo > 0. Solutions with A > 
are inflationary, but without exit, diverging at asymptotic times. Solutions with 
A < 0, independent of a possible inflationary behaviour, asymptotically isotropize 
to a FRW open model. For u 3> 1 the isotropization mechanism happens after the 
(short) inflationary era (see figure [I]), and for to <C 1 it occurs during inflation; some 
figures of the latter case are shown elsewhere [37]. For the radiating case iy = 1/3), 



anisotropic solutions are found for A 7^ and A = 0, as well. Solutions with 
A > are inflationary again, but with the same properties as in the non-radiating 
case: inflation is only present in one scale factor, and there is no exit of it, i.e. the 
solution asymptotically diverges. In the cases with A < inflation takes place, but 
again to must be negative (u < —3/2) to guarantee a sufficient amount of e-folds 
of expansion (N ~ 68). The reason behind this can be found in Eqs. (0) and @. 
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For uj < —3/2 the sign of some source terms (related to the trace T) and kinetic 
terms are reverse and play an inverse role than normally, making possible to have 
transient inflationary solutions. 

It seems that the most interesting cosmological exact solutions in the BD theory 
are those with a BD parameter different from the desired uj > 500, in fact for negative 
u's. For the latter case there have been reported even more solutions than with 
uj > 0, see for instance [^3|, [L6| - However, values of uj < —3/2 correspond to negative 



scalar field densities in the conformally rescaled Einstein frame (g^ w = G<pg^ w ) fl4 |. 
Therefore, such solutions have to be taken with caution. Furthermore, within the 
BD theory alone is not possible to have reheating in the Universe, since there is no 
potential nor any energy transfer to mass terms. However, reheating is necessary 
after inflation, and therefore, a more general theory than BD is in order. In this 
context, we have analysed the case of an induced gravity theory, where after a mild 
inflation period (due to non-minimal coupling) the Bianchi V model isotropizes, 
afterwards inflation due to the potential of the theory follows and reheating takes 
place |33], [57J. Our solutions may be also of physical interest in a scalar-tensor 
theory, perhaps some kind of scalar-tensor with uj(4>), in which the solutions are to 
be valid during some (cosmic) time stage when the physical theory is mimiced by 
a BD effective action, yet with some other value for uj, not restricted to be bigger 
than 500 at the begining of times. 

Finally, we think that much work has to be done in the field of exact solutions 
in the context of modern cosmology. There are many solutions written in scaled 
variables, whose physics is hidden because of its complexity, and because less effort 
has been put on the analysis of its physical meaning. Nice properties as inflation, 
graceful exit, and isotropization has been shown in the present work, and without any 
cosmological constant or function, but for a restricted range of physical parameters. 



References 

[1] Smoot G F et al 1992 Astrophys. J. 396 LI 
Bennett C L et al 1994 Astrophys. J. 436 423 
Wright E I 1994 Astrophys. J. 436 443 
Banday A J et al 1994 Astrophys. J. 436 L99 
Four years of COBE data have been reported by: 
Bennett C L et al 1996 Astrophys. J. 464 LI 
Gorski K M et al 1996 Astrophys. J. 464 Lll 
Hinshaw G 1996 Astrophys. J. 464 L17 

[2] Kolb E W and Turner M S 1990 The Early Universe "Frontiers in Physics" no. 
69 (Addison- Wesley) 



24 



[3] Misner C W 1967 Nature 214 40 

[4] Collins C B and Hawking S W 1973 Astrophys. J. 180 317 

[5] Barrow J D and Sonoda D H 1986 Phys. Rep. 139 1 

[6] Barrow J D 1995 Phys. Rev. D 51 3113 

[7] Gibbons G W and Hawking S W 1977 Phys. Rev. D 15 2752 
Hawking S W and Moss I G 1982 Phys. Lett. B 110 35 

Boucher W and Gibbons G W 1983 The very early universe ed G W Gibbons 
and S W Hawking (Cambridge: Cambridge University Press) 273 

[8] Wald R M 1983 Phys. Rev. D 28 2118 

Steigman G and Turner M S 1983 Phys. Lett. B 128 295 

Futamase T 1984 Phys. Rev. D 29 2783 

Gr0n 1985 Phys. Rev. D 32 2522 

Moss I and Sahni V 1986 Phys. Lett. B 178 159 

Jensen L G and Stein-Schabes J A 1987 Phys. Rev. D 35 1146 

Raychaudhuri A K and Modak B 1988 Class. Quantum Grav. 5 225 

Barrow J D and G6tz G 1989 Phys. Lett. B 231 228 

Kitada Y and Maeda K 1992 Phys. Rev. D 45 1416 

Ibanez J van den Hoogen R J and Coley A A 1995 Phys. Rev. D 51 928 

[9] Anninos P Matzner R A Rothman T and Ryan Jr M P 1991 Phys. Rev. D 43 
3821 

[10] Pimentel L O and Stein-Schabes J A 1989 Phys. Lett. B216 27 

Futamase T Rothman T and Matzner R A 1989 Phys. Rev. D 39 405 

Capozziello S de Ritis R Rubano C and Scudellaro P 1995 Phys. Lett. A201 
145 

Capozziello S and de Ritis R 1995 Phys. Lett. A203 283 
Capozziello S and de Ritis R 1995 Inter. J. Mod. Phys. D5 209 

[11] Kofman L Linde A and Starobinsky A A 1994 Phys. Rev. Lett. 73 3195 

[12] Brans C and Dicke R 1961 Phys. Rev. 124 925 

[13] Kolitch S J and Eardley D M 1995 Ann. Phys., N Y 241 128 
Kolitch S J 1996 Ann. Phys., N Y 246 121 

[14] Holden D J and Wands D 1998 Class. Quant. Grav. 15 3271 

[15] Levin J J 1995 Phys. Rev. D 51 462 

[16] Chauvet P and Cervantes-Cota J L 1995 Phys. Rev. D 52 3416 

[17] Mimoso J P and Wands D 1995 Phys. Rev. D 52 5612 



25 



[18] Will C M 1993 Theory and experiment in gravitational physics (Cambridge: 
Cambridge University Press) chapter 8 

[19] MacCallum M A H 1979 General Relativity. An Einstein Centenary Survey ed 
S W Hawking and W Israel (Cambridge: Cambridge University press) 533 

[20] Ruban V A and Finkelstein A M 1975 GRG 6 601 

[21] Chauvet P Cervantes-Cota J and Nunez- Yepez H N 1991 Proceedings of the 7 
th Latin American Symposium on General Relativity and Gravitation, SILARG 
VII ed J C D'Olivo et. al. (World Scientific) 487 

[22] Wainwright, J and Ellis G F R Dynamical Systems in Cosmology (Cambridge: 
Cambridge University Press). 

[23] Dehnen H and Obregon O 1971 Astrophys. Space Sci. 14 454 
Dehnen H and Obregon O 1972 Astrophys. Space Sci. 15 326 

[24] Lessner G 1974 Astrophys. Space Sci. 30 L5 

[25] Chauvet P and Guzman E 1986 Astrophys. Space Sci. 126 133 

[26] Lorenz-Petzold D 1984 Math. Proc. Camb. Phil. Soc. 96 183 

[27] Lorenz-Petzold D 1984 Astrophys. Space Sci. 98 249 

Ruban V A and Finkelstein A M 1976 Astrofizika 12 371 
Morganstern R E 1971 Phys. Rev. D 4 282 

[28] Barrow J D 1993 Phys. Rev. D 47 5329 

[29] Obregon O and Chauvet P 1978 Astrophys. Space Sci. 56 335 
Chauvet P and Obregon O 1979 Astrophys. Space Sci. 66 515 

[30] Ruban V A 1977 Sov. Phys. JETP 45 629 

[31] Harrison E R 1967 Monthly Notices Roy. Astronom. Soc. 137 69 
Tauber G E 1967 J. Math. Phys. 8 118 
Vajk J P 1969 J. Math. Phys. 10 1145 

[32] Minkowski P 1977 Phys. Lett. B 71 419 
Zee A 1979 Phys. Rev. Lett. 42 417 
Adler S L 1982 Rev. Mod. Phys. 54 729 
Spokoiny B L 1984 Phys. Lett. B 147 39 

Accetta F S Zoller D J and Turner M S 1985 Phys. Rev. D 31 3046 
Pollock M D 1985 Phys. Lett. B 156 301 

Lucchin F Matarrese S and Pollock M D 1986 Phys. Lett. B 167 163 
Pollock M D 1986 Nucl. Phys. B 277 513 

Boutaleb-J. H and Marrakchi A L 1987 Int. J. Mod. Phys. A 2 179 



26 



Accetta F S and Trester J J 1989 Phys. Rev. D 39 2854 

Salopek D S Bond J R and Bardeen J M 1989 Phys. Rev. D 40 1753 

Kaiser D I 1994 Phys. Rev. D 49 6347 

Kaiser D I 1994 Phys. Lett. B 340 23 

Faraoni V 1996 Phys. Rev. D 53 6813 

[33] Fakir R and Unruh W G 1990 Phys. Rev. D 41 1783 
Fakir R and Unruh W G 1990 Phys. Rev. D 41 1792 
Cervantes-Cota J L and Dehnen H 1995 Phys. Rev. D 51 395 
Cervantes-Cota J L and Dehnen H 1995 Nucl. Phys. B 442 391 

[34] Zel'dovich Ya B 1962 Soviet Phys. JETP 14 1143 

Zel'dovich Ya B 1970 Monthly Notices Roy. Astronom. Soc. 160 1 
See, however, Harrison E R 1965 Ap. J. 142 1643 

[35] Obregon O J and Pimentel L O 1978 Gen. Rel. Grav. 9 585 

[36] Turner M S 1983 Phys. Rev. D 28 1243 

[37] Cervantes-Cota J L and Chauvet P A 1999 Phys. Rev. 59 043501 

[38] Levin J J and Freese K 1993 Phys. Rev. D 47 4282 
Levin J J and Freese K 1994 Nucl. Phys. B 421 635 

[39] Magnano G and Sokolowski L M 1994 Phys. Rev. D 50 5039 

Sokolowski L M 1995 Plenary Talk given at the l lth Conference on General 
Relativity and Gravitation, Florence, Italy, re-print |gr-qc/9511073 

[40] Kolb E W Salopek D S and Turner M S 1990 Phys. Rev. D 50 5039 

[41] Fakir R and Unruh W G 1990 Phys. Rev. D 41 1783 
Hwang J 1990 Phys. Rev. D 42 2601 
Guth A H and Jain B 1992 Phys. Rev. D 45 426 
Lyth D H and Steward E D 1992 Phys. Lett. B 274 168 

Starobinsky A A and Yokoyama J 1995 in Proceedings of the fourth workshop 
on gene ral relativity an d gravitation ed K Maeda T Nakamura and K Nakao. 
Reprint |gr-qc/9502002| . 

Mallik S and Chaudhuri D R 1997 Phys. Rev. D 56 625 

[42] Salopek D S Bond J R and Bardeen J M 1989 Phys. Rev. D 40 1753 

[43] Abbott L F and Wise M B 1984 Nucl. Phys. B 244 541 
Lucchin F and Matarrese S 1985 Phys. Rev. D 32 1316 



27 



[44] Lidsey J E et al 1997 Rev. Mod. Phys. 69 373 

Schunck F E and Mielke E W 1994 Phys. Rev. D 50 4794 
Mielke E W and Schunck F E 1995 Phys. Rev. D 52 672 
Mielke E W Obregon O and Macias A 1997 Phys. Lett. B 391 281 
Garcia A Macias A and Mielke E W 1997 Phys. Lett. A 229 32 
Benitez et al 1997 Inter. J. Mod. Phys. A 12 2835 



28 



Figure captions 



Figure 1: The upper graph shows the Hubble parameters as a function of the 
time 7]. As the universe evolves the three Hubble rates tend to their corresponding 
open FRW solution H\. For these plots we have taken v — 0, u — 500, and 
h 2 = r] = m = 1. For these values one has that Hi a > 0, H 3o > 0, and H 2o = 0. 
In the figure below the evolution of the scalar field is shown. For rj ^> 2.4, <p ~ l/r/ 
as in the FRW cosmology. 

Figure 2: The upper graph shows the Hubble parameters as a function of the 
time rj. In these plots we have taken the parameters to be the same as in figure [3], 
except for the anisotropy parameter that we have augmented to be h 2 = 10. Because 
of this, a 2 undergoes a contraction that diminishes as time evolves. At 77 ~ 9 the 
contraction stops and an expansion follows, ai and 03 always expand. The figure 
below shows the scalar field as a function of the time rj. 

Figure 3: The upper graph shows the Hubble parameters as a function of the 
time rj. For these plots we have taken v — 0, — f) — m — 1, and to = 0.1. 
That is, we have changed u to be a small number to show an example of the type 
of induced gravity without potential. The evolution of the scalar field is shown in 
the figure below. 

Figure 4: The upper graph shows the Hubble parameters as a function of the 
time rj for the parameters to = —2, h 2 = c\ = c 2 = -m ]/3 = 1. For these values the 
three initial Hubble functions are positive, Hi Q > 0. The expansion begins super- 
inflationary until some turning point after which the three Hubble rates tend to the 
isotropic solution. In the figure below the evolution of the scalar field is shown. 

Figure 5: The upper graph shows the Hubble parameters as a function of the 
time rj for the parameters uj = —2, h 2 = c\ = 100, c 2 = 10, m 1/3 = —1. For 
these values H\ Q > 0, H 2o > 0, and H 3o < 0. Hi and H 2 begin expanding, and H s 
contracting, until some turning point after which the three Hubble rates tend to the 
isotropic solution. In the figure below the evolution of the scalar field is shown. 

Figure 6: The number of e-foldings (N) of exponential expansion as a function of 
h 2 and c 2 . This plot corresponds to an integration time from rj — to 5. The plane 
N = 68, representing the threshold of successful expansion, is shown for comparison. 
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